Bounds on Spectral Dispersion from Fermi-detected Gamma Ray Bursts 
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Data from four Fermi-detected gamma-ray bursts (GRBs) is used to set limits on spectral dis- 
persion of electromagnetic radiation across the universe. The analysis focuses on photons recorded 
above 1 GeV for Fermi detected GRB 080916C, GRB 090510A, GRB 090902B, and GRB 090926A 
because these high-energy photons yield the tightest bounds on light dispersion. It is shown that 
significant photon bunches in GRB 090510A, possibly classic GRB pulses, are remarkably brief, an 
order of magnitude shorter in duration than any previously claimed temporal feature in this energy 
range. Although conceivably a > 3u fluctuation, when taken at face value, these pulses lead to an 
order of magnitude tightening of prior limits on photon dispersion. Bound of Ac/c < 6.94 x 1CP 21 
is thus obtained. Given generic dispersion relations where the time delay is proportional to the 
photon energy to the first or second power, the most stringent limits on the dispersion strengths 
were ki < 1.61 x 1CT 5 sec Gpc" 1 GeV" 1 and k 2 < 3.57 x 1CT 7 sec Gpc" 1 GeV~ 2 respectively. Such 
limits constrain dispersive effects created, for example, by the spacetime foam of quantum gravity. 
In the context of quantum gravity, our bounds set M\c 2 greater than 525 times the Planck mass, 
suggesting that spacetime is smooth at energies near and slightly above the Planck mass. 

PACS numbers: 98.70.Rz, 11.30.Cp, 98.80.Qc, 14.70.Bh 
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Gamma-ray bursts (GRBs) are the furthest known ex- 
plosions in the universe. Their rapid variability and great 
distances make them useful as probes of light proper- 
ties as well as the intervening space. Were light to have 
fundamentally different speeds at different wavelengths 
(spectral dispersion), distant GRBs might show persis- 
tent energy-dependent arrival patterns [l|. Spacetime 
foam inherent in some formulations of quantum gravity, 
for example, might cause spectral dispersion 0-0] • Other 
properties of light or the universe might also cause differ- 
ent wavelengths to propagate at different speeds [H, 0] . 

GRBs have already been used to limit the cosmologi- 
cal density of compact objects through the non-detection 
of their gravitational lensing [7J. Lag- minimizing al- 
gorithms have been previously designed to search for 
quantum-gravity based dispersion effects Q. Although 
bounds on quantum gravity dispersion in Fermi GRBs 
have been explored previously for two different Fermi 
GRBs 0,0, the present work limits more general pa- 
rameters, considers four Fermi GRBs, considers only 
super-GeV photons, and yields substantially tighter 
bounds. 

Given that two photons of different energies AE are 
emitted at the same place and time, the gap At between 
their arrivals can be quantified as 



At = kr,Dr,E n AE, 



(1) 



where k n is the dispersion strength and D n is a cosmo- 
logical lookback distance that also depends on the nature 
of the photon dispersion [ll| . Specifically, 



Dn = T7- 



(1 + z') n dz' 



H Jo yn^TT^pTfil' 



(2) 



cosmological constant density [11 
flat universe 
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12J in a geometrically 



where H a is present value of Hubble's constant, and VLm 
and tt\ are the present values of the matter density and 



_ For clarity and following theoretical precedents [11|, [14|, 
15| . only three cases will be considered here: n = — 1, 
n = 1 and n — 2. The first case, n — —1, is for a 
universe with no chromatic dispersion. Then, = 
and D_i corresponds with the classic cosmological look- 
back distance [121 ]. In the second case, n = 1, the dis- 
persion delay scales with the energy difference between 
photons, a primary case expected were spacetime to have 
the foaminess inherent in some models of quantum grav- 
ity [l^]. The third case has n = 2, is considered in some 
models of quantum gravity [lij . It will be assumed here 
that dispersion occurs uniformly along the light paths. 

For a group of photons emitted over a source of finite 
size, an upper limit on At might relate primarily to an 
upper limit on source size and not to dispersion proper- 
ties of light. Given limited information, one might not 
be able to disentangle the various contributions to At. 
Surely, though, an observed bound on At would constrain 
the combined processes, thereby limiting the individual 
magnitudes. An exception to this would be if the source 
and universe dispersion effects were of similar magnitudes 
but of opposite sign, a coincidence that is testable with 
a larger data set but here considered unlikely. 

Because the largest energy ranges occur most com- 
monly in the GRBs with the highest energy photons, 
and since these GRBs with many high energy photons 
are rare, GRBs with numerous high energy photons were 
initially sought - to find the finest temporal feature of 
statistical significance. A useful previous search included 
one by Rubtsov et al. 16] of the Fermi LAT photon 
database, although other previous studies also were in- 
fluential 0, [13, flq | ■ Another clue came from a visual 
inspection of Figure 1 of Abdo et al. (2009) [lpj . where a 
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FIG. 1: Time series of photon arrivals for the four Fermi 
GRBs analyzed. Top panel is a closeup of the one second of 
GRB 0901510A, containing the finest temporal features. 

striking clustering of photons above 1 GeV was spotted 
for GRB 090510A. Other reasons for our 1 GeV thresh- 
old include the lower photon background at higher ener- 
gies, and the possibility of extremely brief GRB pulses 
at higher energies. Four candidate GRBs eventually 
emerged: GRB 080916C, GRB 090510A, GRB 090902B, 
and GRB 090926A. "Pass 7" data from these GRBs were 
downloaded from the Fermi web interface at NASA's 
GSFC in February 2012. Only photons within a 95 % 
energy-dependent error radius of the sky position of the 
optical counterpart were considered. This error radius 
was interpolated from Fermi performance data given by 
Ref. 0. 

The bottom four panels of Fig. (1) show time series 
for the arrivals of photons of the four Fermi GRBs. The 
origin t — indicates the time that the GRB triggered 
on Fermi's GLAST Burst Monitor (GBM). On the left, 
at negative times, is 100 seconds of Fermi LAT data that 
occurred before the trigger time while on the right, at 
positive times, is 100 seconds of data that occurred after 
the trigger time. Individual counts are shown as vertical 
line segments. The height of the line segment indicates 
the recorded energy of the photon detected. Inspection 
of Figure 1 shows that the background for stray photons, 
prior to the trigger time, for example, is very low. The 
top panel of Fig. (1) shows a closeup of the one second 
of GRB 0901510A when the bunched photons arrived. 

For a (short) GRB 090510A, consider the hrst 11 pho- 
tons arriving over a AT = 0.1745 seconds. The post- 
trigger arrival times of these photons were 3.702234, 
3.702783, 3.706941, 3.719431, 3.763108, 3.764177, 
3.799190, 3.799318, 3.800096, 3.816729, 3.875767 sec, 
respectively. For comparison, the next five photons, 
photons 12 through 16, arrived at 3.925311, 3.953093, 
4.037660, 4.140611, and 4.152783 sec. The sixth photon 
had the unusually high energy of 30.9 GeV. Of the 11 
photons considered, 6 photons arrived before the tempo- 



ral midpoint and 5 photons arrived later. Notable is the 
closeness in arrival times of three photon groups. These 
groups are defined by the first and second photons, the 
fifth and sixth photons, and photons seven through nine. 
The time between the first and last photons in these 
groups are 0.549 ms, 1.069 ms, and 0.906 ms respectively. 

Is this arrival pattern of remarkably brief doublets sep- 
arated by long pauses significant? Do these 3 brief pulses 
define the finest time scale yet? We argue that such 
"rhythm" is, most likely, not spurious. As shown be- 
low, this group of 11 photons is consistent with a con- 
stant overall arrival rate. Yet, the following simple, al- 
beit crude, analytical argument shows that the odds of 
a uniformly emitting source producing the pattern de- 
scribed above, are below 3 a. This is then confirmed by 
a detailed Monte Carlo simulation. 

For a perfectly random (Poisson) process, the wait- 
ing times (t) between consecutive photon arrivals are ex- 
ponentially distributed and a sum of m such times is 
T-distributed with exponent m (convolution of m expo- 
nential variates) . Given an estimated mean waiting time 
t = 0.1745/10 sec, the probability of waiting t << t is 
t/r. For example, consider t < 1.069 ms a "success" . The 
probability of success is then ps 0.1069/1.750 = 0.0613 
for the 11 photon group. Then the (binomial) proba- 
bility of at least 4 "successes" in 10 trials (10 waiting 
times between the 11 photons) is P(4, 11) = 1/455. If 
one counts the triplet as 3 successes, the odds drop to 
P(5,ll) = 1/5000. These crude estimates bracket the 
result of the 10 9 uniformly random Monte Carlo runs, 
indicating that the chance that 5 photons would trail 
other photons by 1.069 ms or less occurs in only about 
1 in 1190 trials (about 3.34 a). A sceptic might object 
that the mean rate need not be uniform, that both the 
1st and the 11th arrivals ought to be regarded as fixed, 
etc. To that end, we now describe our data analysis as 
well as more elaborate Monte Carlo simulations in more 
detail. 

To determine the briefest yet statistically meaningful 
time interval At in the data, we proceeded as follows. 
Groups of consecutive photon arrival times were consid- 
ered, starting from the three photons arriving closest in 
time, then the four closest photons, and subsequently all 
numbers of GRB-associated photons for 500 seconds fol- 
lowing the trigger. To ensure relatively uniform average 
arrival rates, we chose photon groups with roughly equal 
numbers of photon arrivals before and after the temporal 
midpoint of the group. Formally, a two-bin x 2 statistic 
was computed. Given the single degree of freedom, "flat" 
groups with \ 2 < 1 were considered as statistically con- 
sistent with a flat distribution, and then search for At 
proper ensued, aided by a Monte Carlo simulation as fol- 
lows. 

For each photon in the time series except the last, the 
number of trailing photons, arriving within a time win- 
dow T was counted, for a wide range of Ts. This photon 
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TABLE I: Measured parameters of selected high energy Fermi GRBs. 



GRB 


At 


N 


E (low) 


E (high) 


2 [a][Ref] 


Name 


sec 




GeV (2a) 


GeV (2a) 


(2a) 


080916C 


37.9 


14 


1.32 


10.6 


4.05 [23J 


090510A 


0.00155 


11 


1.58 


24.7 


0.897 [10] 


090902B 


23.9 


33 


1.20 


9.02 


1.82 [24] 


090926A 


3.00 


7 


1.38 


2.15 


2.106 [25] 



TABLE II: Derived and limited parameters of selected high energy Fermi GRBs. 



GRB 




Ac/c 


Dj 


fci 


M lC 2 


D 2 


k 2 


M 2 c 2 


Name 


Gpc 




Gpc 


sec/(Gpc GeV) 


GeV 


Gpc 


sec/ (Gpc GeV 2 ) 


GeV 


080916C 


3.57 


1.03E-16 


16.8 


2.42E-01 


4.25E+17 


48.8 


6.28E-03 


4.96E+09 


090510A 


2.17 


6.94E-21 


4.18 


1.61E-05 


6.41E+21 


6.09 


3.57E-07 


6.57E+11 


090902B 


2.96 


7.84E-17 


8.38 


3.65E-01 


2.82E+17 


16.0 


1.70E-02 


3.01E+09 


090926A 


3.10 


9.40E-18 


9.59 


4.05E-01 


2.54E+17 


19.5 


7.41E-02 


1.44E+09 



count was compared to that expected from a uniformly 
random arrival time distribution. The comparison distri- 
bution typically involved 10 6 trial time-series. To avoid 
spurious bunching, only At delays such that the associ- 
ated number of real photons was found in less than 1 % 
of the equivalent Monte Carlo distributions, were consid- 
ered for further analysis. 

Returning to the GRB 090510A 11 photon group, do 
the 3 brief pulses define the finest time scale of signif- 
icance? Could a variable mean rate, perhaps, produce 
such a pattern? To that end, we assumed that each of 
the 3 photon groups was randomly chosen from a single 
parent pulse form. This pulse form is the generic GRB 
"Norris pulse" shape first suggested by Norri s 12011 for the 
instance found most common by Nemiroff [21] . specifi- 
cally, P = Ae^ t l a " a l t where P is the photon count rate, 
A is the pulse amplitude, t is time during the pulse, and 
a is the time scale of the pulse. To be conservative, we 
will focus on the broadest photon group, the central pair 
separated by 1.069 ms. 

A simple simulation shows that randomly chosen pairs 
of photons from a Norris pulse form have a mean pair sep- 
aration of about 1.20 a. Additionally, in a Norris pulse, 
68.2 % of the photons arrive within a total time window 
of 1.74 a surrounding the pulse peak, here called the 
pulse "width" . Therefore, a parent pulse with width of 
At = (1.74/1.20)1.069 ms = 1.55 ms would yield a mean 
pair separation of 1.069 ms, the longest time between 
first and last photons of the three close photon groups 
of GRB 090510A. Therefore, in subsequent analysis, we 
will use At — 1.55 ms. 

The conservative value of At estimated above for GRB 
090510A is about a factor of ten smaller than even the 
least conservative limit on At listed by Ref. [To[ in row 
5 of Table SI. A primary reason for this is that Ref. [13] 
measured the limiting At essentially as the time differ- 
ence between the start of a sub-MeV spike and a possibly 



associated 0.75 GeV photon. Our analysis differs from 
this earlier analysis of GRB 090510A in that they looked 
at photons over a wide range of energies, whereas we 
looked at only the most energetic photons (>1 GeV) be- 
cause the pulse durations are known to decrease greatly 
as photon energy is increased, so the tightest limits on 
the dispersion delays will come from the highest energy 
photons. Therefore, the small At values presented here 
focus on extremely short doublets prominent at very high 
energies. 

Of the four GRBs considered, only GRB 090510A and 
GRB 0900902B have photons arriving close enough in 
time to eclipse the 0.01 sec previously reported [l(| as the 
smallest At record. We therefore conclude that analyzing 
the other GRBs at most increased the number of trials 
to two, which would decrease the statistical significance 
of the At reported here for GRB 090510A to about 3.14 
a, still above 3 a. 

For GRB 08016C, GRB 090902B, and GRB 090926A, 
none of the photon groups for which the 2-bin y 2 test 
was less than unity showed significant bunching on any 
time scale. On longer time scales, clearly distinct photon 
groups have their At values recorded in Table 1. 

Table 1 lists the measured parameters for the four 
GRBs selected. Column 1 lists the title of the GRB, 
coded with its date of detection. The At values as well 
as the number of photons N on which they are based as 
listed in Columns 2 and 3 respectively. 

Another measured parameter that limits spectral dis- 
persion is AE, the energy between the highest and lowest 
energy photons arriving from the GRB in the At time 
window. Conservative 2 a values of the lowest and high- 
est energy photons - E (low) and E (high) - are given, as- 
suming a 10 % single a energy measurement uncertainty. 
They are listed in Columns 4 and 5 of Table 1 . Values for 
the GRB redshifts were obtained by others from follow- 
up observations of the GRB optical afterglows and the 
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2cT lower limits are listed in Column 6 of Table 1, with 
references. 

The ratio of At and AE has been used to set limits on 
Lorentz invariance previously, where Boggs et al. [22| de- 
rived an upper limit of At/AE of 0.7 sec / GeV for GRB 
021206. For GRB 090510A, Ref. 10] list At/AE < 0.03 
sec / GeV at the 99 % confidence level as their conser- 
vative limit (no least conservative limit is listed). The 
tightest bound from Table 1, however, involving the up- 
per limit on At for GRB 090510A, is At/AE < 6.71 
x 10~ 5 sec / GeV, an improvement of greater than two 
orders of magnitude. 

From the measured parameters listed in Table 1, de- 
rived and limited parameters were computed and listed in 
Table 2. Values of -D-i, Di, and D 2 were computed from 
Eq. (2) under the assumption of a flat concordance cos- 
mology with n,M = 0.3, J7a = 0.7 and a Hubble constant 
H Q of 72 km sec -1 Mpc" 1 , and are listed in Columns 2, 
4, and 7 of Table 2, respectively. 

Given the above data, it is possible to place bounds 
for the difference between the speeds of light at differ- 
ent energies: Ac/c. Assuming Ac results from an in- 
herent property of electromagnetic radiation itself, then 
the lookback distance each photon has traveled is D-\ 
as given by Eq. (2) [13]. Defining lookback time as 
t = D-i/c, the time differential yields Ac/c = cAt/D_±. 
Limits on Ac/c, computed using our strictest upper limit 
on At, are listed in Column 3 of Table 2. 

A previous limit on Ac/c using GRBs was obtained in 
1999 by Schaefer [3], where an analysis of GRB 930229 
yielded Ac/c < 6.3 x 10~ 21 for photons of energies be- 
tween 30 and 200 KeV. A comparable limit for Ac/c < 
6.94 x 10 -21 is derived here from the At listed in Col- 
umn 2 of Table 1 for GRB 090510A for photons of energy 
difference AE > 23.5 GeV. 

Alternatively, it can be assumed that it is the inter- 
vening space that causes differential speed for photons of 
different energies. Following Eq. (2) and approximating 
E ~ AE, it is clear that k n < At/(D n AE n ). In other 
words, were k n greater than this, the universe would have 
separated photons of an energy difference greater than 
AE by more than At. For n = 1 and n — 2, using the 
At limits listed in Column 2 of Table 1, limiting k\ and 

values are listed in Table 2's Columns 5 and 8 respec- 
tively. 

The ki parameter effectively limits dispersion expected 
in some versions of quantum gravity [15j . In particular, 
given that At ~ (AE/Mic 2 )(Di/c) as delineated in Ref. 
[l0(, then Mic 2 = (fcic) -1 . In this parametrization, Mic 2 
is a minimum energy scale of the inherent foaminess of 
spacetime responsible for the dispersion. Note that the 
above data places an upper limit on k\ which translates 
into a lower limit on Mic 2 . Similarly, it is found that 
Af 2 c 2 = (3fc 2 c/2)~ 1 / 2 . The limiting values of Mic 2 and 
M 2 c 2 are listed in Table 2's Columns 6 and 9 respectively. 

Prior to Fermi, GRB published lower limits for 



M 1 /Mpi anck and M 2 /M P i anck were on the order of 0.04 
and 4 x 10~ 12 respectively 14, 22|, where Mpi anc kC 2 = 
1.22 x 10 19 GeV. Using Fermi data for GRB 090510A, 
however, Ref. [1(3] found M\j 'Mpianck > 102, while this 
was relaxed to Mx/Mpi anc k > 1.19 for more conservative 
assumptions. Note that using the most stringent upper 
limit on At listed found here for conservative assumptions 
results in a rather tight bound of Mx/Mpi anc k > 525, 
suggesting that space is smooth even at energies near 
and slightly above the Planck mass. 
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